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Abstract. In this talk I review the status of our ability to extract the CKM matrix elements \Vub\ 
and \Vcb\ from semileptonic decays. I will review both exclusive and inclusive methods and put a 
strong emphasis on how to ensure keeping the extractions model independent. 


INTRODUCTION 

One of the main goals in particle physics over the next several years is to test precisely 
the flavor sector of the standard model (SM). In the SM, all weak interactions are 
determined by 12 parameters, namely the Fermi coupling constant, the weak mixing 
angle Bw, the six quark masses and four parameters in the CKM matrix. B physics plays 
an integral role in overconstraining the four parameters in the CKM matrix, which leads 
to a stringent test of the SM. In this talk I will give an overview of the progress that 
has been made on using semileptonic B decays to determine two of these parameters, 
the magnitude of Vut and Vet- Semileptonic decays provide an ideal way to measure 
the magnitude of CKM matrix elements, since the strong interaction effects are greatly 
simplified by the presence of the two leptons in the final state. 

I will review how to measure Vub and Vcb from both inclusive and exclusive decays. 
I will present the current status of extracting these CKM parameters from semileptonic 
decays, and point out possible improvements in the coming years. 


DETERMINATION OF Kb 
Exclusive Decays 

One of the main applications of heavy quark effective theory (HQET) [1] is that Vcb 
can be extracted in a model independent way using exclusive b^c semileptonic decays. 
Heavy quark spin and flavor symmetry predicts that in the infinite mass limit all heavy to 
heavy form factors are determined by a single nonperturbative function, the Isgur-Wise 
function. Moreover, this function is normalized to unity at zero recoil, which allows for a 
model independent prediction of semileptonic exclusive decays at that kinematic point. 
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To be more precise, consider the decay B D*iv, and define the four-velocities of 
the B and the D* mesons 


^h~ ^ 1 




Pd* 


mg “ mg)* 

The decay rate as a function of the recoil w = ■ Vc is given by 


dTg^D* 
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If {mg, mo, w) I Vcb \ (w), 


where — 1 is a phase space suppression factor, (w) is a form factor and 

N S/, n?/, 7 n 4w l—2wr + r^\ 

f(,n„ w) = (1 -.) (1 + w ) ( 1 + (i _,)2 j . 

where r = mg* j mg. The form factor is related to the Isgur Wise function by 

Fd*{w) =riA^(w), 
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where is a correction factor including both perturbative corrections, given by an 
expansion in <Xs{m},), and nonperturbative corrections, determined by an expansion in 
Almb^c- At zero recoil the Isgur-Wise function is normalized, ^(1) = 1, and all 0{A/m) 
corrections in vanish at that kinematic point [2]. The perturbative contributions to 
have been calculated to two loop order [3] and the A^/corrections can be estimated 
using phenomenological models [4]. This leads to [5] 

Fz)*(l) =riA = 0.913 ±0.042, (5) 

where the uncertainty is mostly due to the unknown A^/m^ and higher corrections. 
Recently, this form factor has also been calculated on the lattice, and the result has 
identical central value and comparable uncertainties^ [6] 

FLf(l)= 0.913^24+17. (6) 

Several experiments have measured this differential decay rate and for illustration 
we show in Fig. 1 a plot of the BELLE collaboration together with a plot comparing 
the different measurements. Since the rate vanishes at zero recoil due to the phase 
space suppression factor, the differential rate has to be extrapolated to zero recoil. This 
extrapolation is guided by a model independent relationship between the curvature and 
the slope of ^£>*(1) [7]. Combining the experimental measurement with the theoretical 
calculation of Fo*(l) leads to the values presented at this conference [9, 10, 11] 


t/BELLE 
^cb — 

[39.7 ± 1.6{stat) ±2.0(5y5) ± l.9{theor)] x 10 ^ 


vCLEO 
^cb — 

[46.2± l.4{stat) ±2.0(5y5) ±2.1{theor)] x 10 ^^ 

(7) 

= 

[40.5 ±0.5 ±2.0] X 10^2 

(8) 



(9) 


7 This result does not include the QED corrections of -1-0.007, which is included in (5). 




FIGURE 1. Left; Extrapolation of \Vcb\FD* (w) to the zero recoil point by Belle. Right: comparison of 
results from various collaborations. 


I would like to make a few eomments on how the aeeuraey in this measurement might 
be improved. As ean be seen from Fig. 1, there is a large eorrelation between the slope of 
the physieal form faetor and its value at zero reeoil. Thus, information about this slope 
ean reduee the errors on \Vcb\ eonsiderably. There have reeently been suggestions that 
QCD sumrules give eonstraints on the slope of the Isgur-Wise funetion [12], although 
one has to be eareful in relating these to the slopes of the physieal form faetor. Sinee the 
slopes of the B D and the B ^ D* form faetor are related to one another in a ealeulable 
way, one should use the information from both of these deeays at the same time (note 
that the experimental uneertainty in the slopes in these two deeays is roughly the same 
[9]). In summary, performing a new eomplete high statisties analysis, ineluding both the 
D and D*, as well as the form faetor relations Ri and R 2 , should lower the uneertainties 
below the eurrent level of 5%. 


Inclusive Decays 


Inelusive semileptonie B deeays ean be ealeulated using an operator produet expan¬ 
sion (OPE). This leads to a simultaneous expansion in powers of the strong eoupling 
eonstant CXj(m^) and inverse powers of the heavy b quark mass. At leading order in this 
expansion this reproduees the parton model result 


ro = (1 - 8p + 8p3 - p" - 12p"logp) . (10) 

where p = m^/m^, and nonperturbative eorreetions are suppressed by at least two powers 
of nib. The resulting expression for the total rate of the semileptonie B —0 Xciv ean be 
written as 
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where the three eoeffieients A, B,C depend on the quark masses . The perturbative 
eorreetions are known up to order a^(3o. There are no nonperturbative eorreetions at 



















order Ajnit,, and at order A^jm^ they are given in terms of the two matrix elements 
Xi and X 2 , with the dependenee on these matrix elements contained in the coefficient 
C = C(?ii, X 2 ). Since the value of X 2 can be obtained from the B — B* mass splitting, and 
the charm quark mass can be related to known meson masses, the b quark mass and Xi, 
the only unknowns in extracting \ Vcb\ from the inclusive rate are the precise value of the 
b quark mass nib and the matrix element Xi. 

In my opinion the best way to extract nib and Xi is to use the semileptonic decay data 
itself. Various differential decay spectra can be measured, not only the total decay rate. 
Several observables have already been constructed out of these spectra which are sensi¬ 
tive to nib and Xi [13]. The goal should be to measure and calculate as many inclusive 
observables as possible. All these observables are given in terms of nib and Xi, so a fit 
to these observables will determine these parameters with theoretical uncertainties given 
by A^ jm^ terms in the OPE [14]. Recently, the CLEO collaboration performed such an 
analysis in which they used the first moment of the hadronic invariant mass spectrum in 
B Xc£v and the first moment of the photon energy in the rare decay B Xsj to extract 
nib and Xi and used these values to determine V^b [10]. Defining A = nig — nib +... they 
find 


A = 0.35±0.07±0.10GeV, Xi =-0.238 ±0.071 ±0.078 GeV^ (12) 
and with that 

\V,b\ = [4.04±0.09±0.05±0.08] x 10~^ (13) 

The errors are (in order) from the measurement of the total decay rate, from the uncer¬ 
tainties in A, Xi, and from the higher order terms in both the perturbative and nonpertur- 
bative expansion. 

Using A as the mass parameter in the semileptonic decay rate has the disadvantage 
that the parameter A is not an infrared safe quantity. This means that the value of this 
parameter is ambiguous up to terms of order Aqcd [15]. This ambiguity cancels in any 
physical observable by a similar ambiguity in the perturbative expansion [16]. Eor this 
cancellation to occur, however, it is essential that one works consistently to a particular 
order in perturbation theory. Since this is an artifact of the pole mass, it would be 
convenient to repeat the analysis with an infrared safe definition of the b quark mass 
[17]. The numerical value of such a short distance mass can then be used safely in other 
processes as well. 

Einally, I want to comment on an uncertainty not included in (13), namely the uncer¬ 
tainty from duality violations, which are present for any inclusive observable in heavy 
quark decays. While there have been several estimates of these uncertainties [18], I again 
would like to advocate using the semileptonic decay data itself to test for duality viola¬ 
tions. Since many semileptonic observables can be calculated to the same order as (11) 
and all depend on the same two parameters nib and Xi, one can check for duality vi¬ 
olations by comparing the theoretical predictions with experimental results. This is an 
important analysis which still needs to be performed. 



DETERMINATION OF Kb 


Exclusive Decays 


Measuring \ Vub \ from exclusive decays of B mesons is considerably more difficult than 
the determination of \Vcb\ discussed earlier. This is because the heavy quark spin-flavor 
symmetry, which gives rise to the simplifications of heavy to heavy form factors, does 
not apply to the heavy to light decays mediated by a Z? ^ m transition. Thus, to extract 
114^1 from exclusive decays one has to find other ways to calculate these form factors. 

As an example, consider the exclusive decay B Ttfv. The required matrix element 
can be written in terms of two form factors 


{n\ufPL\B) = U{q^) 


{PB+Pn) 






9 H : 
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(14) 


and a similar relation is true for B —> pfv. The form factor /o(^^) vanishes for zero 
lepton mass, but the form factor f+{q^) is required to describe the decay rate. There are 
various suggestions in the literature how to extract these form factors. For B pfv 
the relevant form factors can be related to those in D — K*(m using heavy quark 
and chiral symmetries [19]. Corrections are first order in both heavy quark and chiral 
symmetry breaking individually and are at the 30% level. This uncertainty could be 
reduced considerably if information from D pfv and B became available, 

since then a double ratio can be constructed which has corrections which simultaneously 
violate heavy quark and chiral symmetry [20]. The resulting uncertainties should be 
below the 10% level. 

Ultimately, lattice calculations will determine the required heavy to light form factors 
in a model independent way. In the last year several quenched calculation of the form 
factor f+{q^) have been performed. Since lattice calculations can only be performed 
for slow pions, the lattice results are only available for small pion energies E,! <,1 GeV. 
A result from a recent lattice calculation of this form factor [21] is shown in Fig. 2. 
Experimentally, the situation is reversed with the efficiency to measure slow pions being 
small. In order to use lattice calculations to determine \Vub\ from exclusive decays, 
however, it is crucial to increase the kinematic overlap of the lattice calculation of the 
form factors with the experimental measurements of the differential decay rates. 


Inclusive Decays 

The inclusive decay rate B —> X^iv is directly proportional to and can be 

calculated reliably and with small uncertainties using the operator product expansion. 
Unfortunately, the ~100 times background from B Xc£v makes the measurement of 
the totally inclusive rate an almost impossible task. Several cuts have been proposed 
in order to reject the b ^ c background, however care has to be taken to ensure that 
the decay rate in the restricted region of phase space can still be predicted reliably 
theoretically. The cut which is easiest to implement experimentally is on the energy of 
the charged lepton, requiring > (m| — m^,)/2mg. Unfortunately, this cut restricts the 




FIGURE 2. Lattice calculation of the form factor /+ as a function of the pion energy. 


remaining region of phase space too much for the OPE to still be valid. Instead, a twist 
expansion has to be performed [22], and at leading order the decay rate is determined 
by the light cone distribution function of the B meson, with subleading twist corrections 
suppressed by powers of A/m^, [23]. This distribution function is the matrix element 
of a nonlocal operator and can not be calculated perturbatively. However, since it is 
a property of the B meson itself, and thus independent of the decay process, it can 
be extracted from a different process and then used in the inclusive b —> decay. 

This is similar to parton distribution functions, which can be measured for example 
in ep collisions and then be used in other processes. The best way to measure this 
distribution function is in the decay B Xsj. At leading order in a^, the shape of the 
photon energy spectrum is entirely due to this light cone distribution function, and the 
perturbative corrections are known to order [24]. The photon spectrum has recently 
been measured to good accuracy by the CLEO collaboration [10], which constituted the 
first determination of this structure function. 

Unfortunately, a cut on the lepton energy has the disadvantage that only ~ 10% of 
the events pass this cut. Thus, the decay is far from being fully inclusive and one has to 
worry about duality violations becoming large. It is therefore advantageous to find ways 
to suppress the background which allows more b ^ u events to survive. The optimal 
such cut is on the hadronic invariant mass nix < [25]. This cut is optimal in the sense 

that any other cut is a subset of it, and it has been estimated that ~ 80% of the b ^ u 
events survive. The same structure function which describes the lepton energy endpoint 
is needed to describe the region of low invariant mass. In [26] it was shown how to 
eliminate the structure function entirely and measure \ Vub\ model independently by using 
semileptonic data in combination with the photon energy spectrum in 5 —> X^y. The 
main uncertainties to this method of measuring \ Vub\ are given by unknown subleading 
structure functions, which are parametrically suppressed by Aqcd/ mb [23]. It is hard to 
quantify the exact size of these uncertainties, and I would conservatively estimate them 
to be 10-20%. 

It was shown in [27] that using a cut on the leptonic invariant mass ~ {p£ + pv)^ 
allows to measure \Vub\'^ without requiring knowledge of the structure function of the 
B meson. The number of events surviving such a cut on can be calculated using the 
usual local OPE and depending on the exact value of the cut chosen, the fraction of 
events surviving the cut is 10-20%, with uncertainties on ranging from 15% for 
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FIGURE 3. The dalitz plot in the /sh and /Ei plane. The grey blob in the lower right (left) hand 
corner of the phase space indicates where the local OPE breaks down. 


^cut = ~ 11.6 GeV^ to 25% for = 14GeV^. The advantage of not having to 

use information from B Xsj and of the absenee of power eorreetions at order Aqcd / i^b 
is offset by the fact that this cut again eliminates a huge fraction of the events and the 
worries about duality violations reappear. 

To improve the situation, one has to understand why the cut allows to calculate 
using the local OPE, while the other cuts require the twist expansion. The twist expan¬ 
sion is relevant if the remaining phase space is dominated by hadronic states with energy 
much larger than their invariant mass. This is the case for both a cut on Eg and a cut on 
mx- The cut on q^, however, restricts the hadron energy to be Ex < mg — \/^ < mo 
and this cut therefore eliminates this dangerous region of phase space. The situation is 
illustrated in Fig. 3. 

In [28] a new strategy to measure \Vub\ was proposed which combines the advantages 
of the nix cut with those of the q^ cut. The idea is to use a cut on the hadronic invariant 
mass to reject the b ^ c background, and simultaneously use a lower cut on q^ to avoid 
the twist region. The pure mx and q^ cuts are contained in this approach as the limits 
^cut = 0 or <?cut = {^B — respectively. Thus, varying the iri the presence of a 

cut on mx allows to smoothly interpolate between these two cases. Using these results 
the strategy to extract \ Vub\ in a model independent way is: 

• make the cut on mx as large as possible, keeping the background from b ^ c small. 

• for a given cut on mx, reduce the q^ such as to minimize the overall uncertainty. 

To illustrate this procedure I will present graphically the results of [28]. In Fig. 4 the 
fraction of events surviving such combined cuts are shown for three values of in 
combination with the effect of the structure function. As advertised, lowering the cut 
on q^ increases the fraction of events, while the effect of the structure function grows. 
Depending on how well the structure function will be known in the future determines 
how much the cut on q^ can be lowered. Due to the fact that the structure function 
can only be determined up to corrections scaling as Aqcd/"*Zi, it will always have an 
uncertainty at the 10-20% level, even with a perfect measurement of the B X^yphoton 
spectrum. Thus, for a determination of \Vub\ with uncertainties below 10% one should 
ensure that the difference between the prediction with and without the structure function 
is less than ~ 30%. 

As in any inclusive observable, the remaining uncertainties are due to three sources: 



























FIGURE 4. Left; as a function of for mcut = 1.86GeV (solid line), 1.7GeV (short 

dashed line) and 1.5GeV (long dashed line). Right; The effect of the structure function. 


unknown perturbative eorreetions, uneertainties due to the parameters and Xi and un- 
eertainties due to unknown matrix elements of loeal operators at 0(A^/m^) in the OPE. 
To eompare various methods to determine \Vub\ it is erueial to have a eonsistent seheme 
to determine the theoretieal uncertainties. In [28] such an error analysis was performed 
for the combined q^, mx cut method, and in Table. 1 I present the result for various dif¬ 
ferent combinations of cuts. From this table it is obvious that even without any knowl- 


TABLE 1. Fraction of events for several different choices of along with the uncertain¬ 

ties on \Vub\ - The two last lines correspond to pure q^ cuts, mcut = nig — and are included for 
comparison. Astmct gives the fractional effect of the structure function f{k+) in a simple model; we 
do not include an uncertainty on this in our error estimate. The overall uncertainty A|V„;,| is obtained 
by combining the other uncertainties in quadrature. The two values correspond to = ±80MeV 
and ±30MeV. 


Cuts on [q^^ m\) 


Fract 


AstructTui 


Apert Vy/? 


^rrih^ub 


±80/30MeV 


iVub 


AVub 


Combined cuts 


6GeV^ 1.86 GeV 

46% 

-2% 

4% 

7%/2.5% 

3% 

8%/5% 

8GeV2,1.7GeV 

33% 

-3% 

6% 

8%/3% 

4% 

9%/6% 

llGeV2,1.5GeV 

18% 

-4% 

13% 

9%/3.5% 

8% 

14%/11% 


Pure q^ cuts 

{mB-mof^mD || 0.17% | — | 15% 19%/7% 18% | 315%/12% 


edge of the structure function, and using a realistic cut of = 1.7 GeV, \Vub\ can still 
be measured with uncertainties below the 10% level by choosing 8GeV^. With 

many different methods to measure \Vub\, the question arises of which measurement to 
ultimately use. In my opinion, this question will most likely be answered once a full 
analysis is performed and the theoretical errors are investigated in conjunction with the 
experimental uncertainties. Depending on the details of the experimental efficiencies, 
some combination (most likely non-linear) of the three cuts on Ei, nix and q^ will prove 
to be the most effective way to minimize the overall uncertainty. The cut on the hadronic 
invariant mass mx will most likely be the best way to eliminate the b ^ c background, 
the cut on the leptonic invariant mass q^ will limit the effect of the nonperturbative struc¬ 
ture function and the cut on the lepton energy Ei can help with experimental efficiencies. 
A careful analysis, which includes a realistic structure function and conservative errors 
should be performed to find the ideal combination of cuts. 




















CONCLUSIONS 


In this talk I reviewed the current status of determining the magnitude of the CKM matrix 
elements \Vub\ and \Vcb\ from semileptonic B meson decays. For exclusive decays I 
reviewed theory and experiment and pointed towards a few improvements possible in the 
coming years. For inclusive decays, I reviewed how their decay spectra can be calculated 
(up to duality violations) using an operator product expansion, which at leading order 
reproduces the parton model results. Nonperturbative corrections are parameterized by 
matrix elements of local operators suppressed by powers of nib. There are no corrections 
at 0{Almb), and at second order in the inverse heavy quark mass there are two matrix 
elements, commonly labeled by Xi and X 2 . While the value of X 2 can be obtained from 
the B — B* mass splitting, the value of Xi is currently not known very well. Thus, in order 
to calculate the inclusive semileptonic observables precisely, we need to obtain values 
for the b quark mass, which enters the calculation, as well as the for the parameter Xi. 

For B Xciv, these parameters give rise to the dominant theoretical uncertainty in the 
determination of \Vcb\ and their precise determination is therefore crucial to lower this 
uncertainty. I have argued here that these parameters should be extracted from decay 
spectra itself, and then used to extract \Vcb\. I also put forward the idea that duality 
violations can be checked using semileptonic decay spectra. 

To measure \Vub\ from the inclusive decay B —one has to deal with the large 
background from b-^ c transitions. Imposing kinematic cuts to suppress this background 
tends to destroy the convergence of the OPE and for both a cut on the lepton energy and 
a cut on the hadronic invariant mass an incalculable structure function is required, with 
corrections suppressed by A/m/,. This function can be extracted from the photon energy 
spectrum in 5 X^y, thus a precise measurement of this spectrum is very important. A 
cut on the leptonic invariant mass can be used to lower the effect of this structure 
function and is therefore crucial to minimize theoretical uncertainties. Ultimately, a 
dedicated analysis which combines the experimental and theoretical uncertainties will 
decide which combination of cuts will allow for the most precise determination of \Vub\. 
I anticipate that such a study will eventually allow a model independent measurement of 
\Vub\ with uncertainties at the 5% level. 
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